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ABSTRACT 

In this work, we investigate the thermophysical properties of OSIRIS-REx target as¬ 
teroid (I0I955) Bennu (hereafter, Bennu), where thermal inertia plays an important 
role in understanding the nature of the asteroid’s surface, and will definitely pro¬ 
vide substantial information for the sampling return mission. Using a th ermophysical 
mode l incorporating the recently updated 3D radar-derived shape model (jNolan et al.l 
l2013f) and mid-infrared observations of Spitzer-PUI, Spitzer- IRAC, Herschel/PACS 
and ESO VLT/VISIR (|M-uller et al.ll20l3 lEmerv et al.l[2M3l . we derive the surface 
thermophysical properties of Bennu. The asteroid has an effective diameter of 5 IOI 4 Q 
m, a geometry albedo of O.OdTto^ooiiJ ^ roughness fraction of 0.04to'o4j and thermal 
inertia of 240^60° Jm ^s ^ for a best-fit solution at Icr level. The best-estimate 

thermal inertia indicates that fine-grained regolith may cover a large area of Bennu’s 
surface, with a grain size that may range from 1.3 to 31 mm, and our outcome further 
supports that Bennu would be a suitable target for the OSIRIS-REx mission to return 
samples from the asteroid to Earth. 

Key words: radiation mechanisms: thermal - minor planets, asteroids: individual: 
(101955)Bennu - infrared: general 


1 INTRODUCTION 

(101955) Bennu (1999 RQ36) is an Apollo-type near-Earth 
asteroid (NEA) from a dynamical viewpoint. As is well 
known, it is widely believed that Bennu could be a potential 
Earth impactor with a rel atively high impact probability of 
appr oximately 3.7 x 10“"^ jMilani et ahllioo^ : ICheslev et al.l 
l2014l 'l. Nowadays, Bennu is recognized as one of the poten¬ 
tially hazardous asteroids (PHA) and its orbit makes it es¬ 
pecially accessible for spacecraft. Therefore, NASA selected 
Bennu as an ideal target for the OSIRIS-REx sample re¬ 
turn mission llLauretta et al.l[2012l l. which will be launched 
in 2016. 

On the other hand, based on a linear, featureless spec¬ 
trum from 0.5 to 2.5 /rm, Bennu is categorized as a B- 
type asteroid, and the related spectral analogue is consid¬ 
ered to be CM chondrite meteorites (!Bus fc Binzell I 2 OO 2 I : 
IClark et al.ll20lih . B-type asteroids are usually thought to 
be primitive and volatile-rich, and most of them are believed 
to originate from the middle and outer main belt. Primitive 
asteroids can offer key clues for us to understand the process 
of planetary formation, the environment of early solar neb¬ 
ulae, and the emergence of life forms on Earth. From spec- 
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troscopic analysis of B-type asteroids, the surface composi¬ 
tion appears to be anhydrous silicates, h ydrated silicates, or¬ 
ganic polymers, magnetite, and sulfid es ( LMSo^^et_^ 198^; 
IClark et al.l l20ld : IZiffer etliD I 2 OIII : Ide Leon et al.l 2012l l. 
whereas those in the outer main belt may also support 
the H 2 0-ice stuff on the subsurface llRivkin fc Emer(1l20ld : 
ICampins et al.ir20ld l. Both spectral and dynamical investi¬ 
gations suggest that Bennu may be a liber ated member of 
the Polana family in th e inner main belt (ICampins et al.l 
l20ld: iBottke et allboisll . 


According to radar observations acquired from Gold- 
stone and Arecibo at its two close st Earth-app r oache s in 
1999 and 2005, the latest analysis of iNolan et al.l (l2013ll re¬ 
solved a nearly spherical shape for Bennu, where its effective 
diameter is 492±20 m, the rotation period is 4.2976±0.0002 
h and spin axis is /3 = —88.0°, A = 45.0°, indicating a 
retrograde rotation. The rotation period is slow enough for 
the nearly spherical asteroid to allow potential sampling of 
the surface regolith. It seems that the rotation period has 
not been greatly spun up by tidal or radiati on forces (e.g. 
YORP effect). However, the recent work of ICheslev et al.l 
(l2014h shows that the semimajor axis of Bennu has a mean 
drift rate da/dt = (—19.0 ± 0.1) x 10““* au • Myr“^ because 
of the Yarkovsky e ffect. On the basis of this estimation, 
ICheslev et al.l (120141 ') predicted numerous potential impacts 
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for Bennu in the years from 2175 to 2196. It is very neces¬ 
sary to continue to observe and study the orbital evolution 
of Bennu. Here we aim to investigate the thermal inertia 
of Bennu, because of a close relationship between thermal 
inertia and Yarkovsky and YORP effects. 

The thermal inertia of an asteroid may be evalu¬ 
ated by fitting mid-infrared observations with a ther- 
mophysical model to reproduce mid-IR emission curves. 
iMuller et al.l ll2012l l estimated Bennu’s thermal inertia to be 
~ 650 Jm“^s“'’ ®K“^ with thermophysical model (TPM), 
based on observations from Herschel/ PACS, ESO-VISIR , 
Spitzer-IRS and Spitzer-PUI. Recently. lEmerv et al.l ll2014l l 
showed that a new result of Bennu’s thermal inertia is 
evaluated to be 310 ± 70 Jm~^s“° ®K“'^, using Spitzer-IRS 
spectra and a mult i-ban d thermal lig h tcurv e. The work of 
I Emery et al.l (120141 ') and iMtiller et al. ll2012ll differs in two 
main aspects: first, the former a dopted a 3D radar-derived 
shape model dNolan et al.ll20f3l in their modeling process, 
whereas the latter utilized a simple spherical shape model in 
the investigation; second, the observational data they used 
were slightly diffe rent, in that IRAC a nd IRS peak-up data 
were included in [ Emery et al.l ll2014l l. but not utilized in 
iMuller et al.l ll2012 b 


In the present work, we adopt ind ependently de¬ 
velop ed thermophysical simulation codes iIYu. Ji &: Wand 
12014 ) based on the fram ework of the Advance d Thermal 
Physical Model (ATPM) dRozitis fc GreerJl201ll) . to inves¬ 
tigate the surface thermophysical characteristics of Bennu. 
In our modelling process, we utilize the rada r-derived shape 
model of asteroid Bennu dNolan et al. 2013|) rathe r than a 
spherical approximation shape ( Mitller et al.l[201^ . In ad¬ 
dition, the mid-infrared data come from four groups of ob¬ 
servations, at various phase angles, by Spitzer- PUI, Spitzer- 
IRAC Herschel/PACS an d ESO VLT/VISIR dMllller et all 
I2OI2I : lEmerv et al.l l2014l ). With the rotational averaged 
ATPM results fitted to all observations, we provide con¬ 
straints for thermal inertia and roughness fraction simulta¬ 
neously at la level. The derived val ue of thermal inerti a 
for Bennu is slightly lower than that of lEmerv et al.l d2014h . 
which may also provide important evidence of the existence 
of fine-grained regolith on Bennu’s surface. Moreover, on 
the basis of the derived thermal inertia, we estimate the 
grain size of the regolith acco r ding t o a thermal conductivity 
model of iGundlach fc BlumI d2013l ). Finally, we summarize 
the major outcomes and present a brief discussion. 


2 THERMOPHYSICAL MODELLING 

2.1 Mid-infrared Observations 

As described above, in our fitting procedure we simply 
adopt the observations from Spitzer- PUI, Spitzer-IRAC , 
Herschel/PACS an d ESO VLT/VISIR (iMwller et al.ll2012l: 
lEmerv et al.l |2014| ). as there are no published data sets 
from Spitzer-IRS. However, the tests of synthetic data from 
Spitzer-IRS were generated, and included for fitting to ex¬ 
amine the difference, but we do not find any significant vari¬ 
ation in the results. 

If the mid-infrared data at various observational phase 
angles, especially at low phase angles, are available, a rela¬ 
tively more reliable thermal inertia may be derived on the 


basis of ID thermal models. However, currently all published 
mid-infrared data of Bennu were observed at large phase 
angles (> 60°), e.g., those from Spitzer-IRS, Spitzer-PUI, 
Spitzer-IRAC, Herschel/PACS and ESO VLT/VISIR. Al¬ 
though the observations at high phase angles constitute a 
disadvantageous condition to constrain thermal inertia and 
surface roughness for the asteroid, we can still use these 
observations, for a combination fitting to observations per¬ 
formed at several different phase angles may remove the de¬ 
generacy of thermal inertia and surface roughness in the 
modelling process, which could somewhat offset the disad¬ 
vantage of the lack of low phase angle observations. From 
this viewpoint, the absence of Spitzer-IRS observations may 
not significantly influence our results, which accord well with 
the aforementioned tests with synthetic data. Hence, we tab¬ 
ulate all data used in the fitting in Table [T] 


2.2 Advanced thermophysical model 

In ATPM dRozitis fc Creeiill201ll : IVu. Ji fc Wandlioidl ). an 
asteroid is considered to be a polyhedron composed of N tri¬ 
angle facets. For each facet, the conservation of energy leads 
to an instant thermal equilibrium between sunlight, thermal 
emission, thermal diffusion, multiple-scattered sunlight and 
thermal-radiated fluxes from other facets. If each facet is 
small enough, the thermal diffusion on the asteroid can be 
approximatively described as one-dimensional (ID) heat dif¬ 
fusion. Hence, the temperature Ti of each facet varies with 
time as the asteroid rotates. In this process, Ti can be en¬ 
larged by multiple-scattered sunlight and thermal-radiated 
fluxes from other facets, which explains well the so-called 
thermal infrared beaming effect. When the entire asteroid 
comes into the final thermal equilibrium state, T will change 
periodically following the rotation of the asteroid. Therefore, 
we can build numerical codes to simulate T at any rotation 
phase for the asteroid. For a given observation epoch, ATPM 
can reproduce a theoretical profile to each observation flux 
as: 

N 

E^odei(A) = Y, Ti) , (1) 

i=l 


where e is the emissivity, f{i) is the view factor of facet i to 
the telescope and B{X,Ti) is the Plank function: 


B{X,Ti) 


2-KhY 1 

exp(^-l) 


( 2 ) 


Thus the calculated Emodei can be compared with the ther¬ 
mal infrared fluxes summarized in Table [T] in the fitting pro¬ 
cess. 


2.3 Fitting Procedure 

In order to derive the thermophysical nature of Bennu via 
the ATPM procedure, we need to know several physical pa¬ 
rameters - the shape model, surface roughness, geometric 
albedo, thermal inertia, thermal conductivity and thermal 
emissivity - which are used as the initial parameters in the 
calculations. In addition, the heliocentric distance, geocen¬ 
tric distance and phase angle of Bennu are well determined 
because the asteroid’s orbit has been accurately measured 
by optical observations. On the other hand, we employ the 
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Table 1. Observational data used in this work 


ifrom [Mriller et alJ ll2012h and lEmerv et ahl |2o3)) 


UT 

16.0 (fim) 

Flux (mjy) 

22.0 

'^helio 

(AU) 

^obs 

(AU) 

a 

n 

Observatory 

Instrument 

2007-05-03 00:00 

12.49±0.28 


11.85±0.41 

1.12380 

0.50564 

-63.52 

Spitzer-PUI 

2007-05-03 02:32 

12.70±0.28 


11.94±0.42 

1.12444 

0.50650 

-63.48 

Spitzer-PUI 

2007-05-03 03:29 

13.38±0.28 


12.40±0.42 

1.12864 

0.51222 

-63.47 

Spitzer-PUI 

2007-05-03 04:40 

12.69±0.29 


12.09±0.40 

1.12516 

0.50748 

-63.46 

Spitzer-PUI 

2007-05-03 08:10 

12.87±0.27 


12.05±0.41 

1.12816 

0.51155 

-63.41 

Spitzer-PUI 

2007-05-03 08:33 

12.47±0.28 


11.83±0.42 

1.12528 

0.50762 

-63.40 

Spitzer-PUI 

2007-05-03 09:55 

13.00±0.28 


12.23±0.40 

1.12414 

0.50611 

-63.38 

Spitzer-PUI 

2007-05-03 10:44 

13.00±0.28 


12.28±0.43 

1.12657 

0.50938 

-63.25 

Spitzer-PUI 

2007-05-03 20:08 

13.08±0.30 


12.18±0.44 

1.12428 

0.50629 

-63.25 

Spitzer-PUI 

2007-05-04 07:40 

12.85±0.29 


12.17±0.42 

1.12492 

0.50716 

-63.10 

Spitzer-PUI 

2007-05-04 11:11 

12.97±0.28 


12.62±0.42 

1.12497 

0.50722 

-63.05 

Spitzer-PUI 


UT 


3.6 (pm) 

Flux (mJy) 

4.5 (pm) 5.8 (pm) 

8.0 (fim) 

^helio 

(AU) 

^obs 

(AU) 

OL 

n 

Observatory 

Instrument 

2007-05-08 

16:24 

0.0545±0.0033 

0.324±0.009 

1.396±0.043 

4.968±0.044 

1.14245 

0.53224 

-61.78 

Spitzer-IRAC 

2007-05-08 

16:46 

0.0482±0.0042 

0.327±0.012 

1.396±0.043 

4.983±0.067 

1.14250 

0.53232 

-61.78 

Spitzer-IRAC 

2007-05-08 

17:14 

0.0481±0.0029 

0.332±0.006 

1.495±0.035 

5.063±0.074 

1.14256 

0.53241 

-61.77 

Spitzer-IRAC 

2007-05-08 

17:39 

0.0488±0.0032 

0.334±0.008 

1.397±0.030 

4.970±0.060 

1.14261 

0.53250 

-61.77 

Spitzer-IRAC 

2007-05-08 

18:06 

0.0494±0.0032 

0.331±0.009 

1.382±0.028 

4.876±0.066 

1.14268 

0.53259 

-61.76 

Spitzer-IRAC 

2007-05-08 

18:27 

0.0450±0.0027 

0.335±0.013 

1.340±0.042 

4.985±0.047 

1.14272 

0.53266 

-61.76 

Spitzer-IRAC 

2007-05-08 

18:59 

0.0425±0.0022 

0.308±0.009 

1.353±0.049 

4.664±0.049 

1.14280 

0.53278 

-61.75 

Spitzer-IRAC 

2007-05-08 

19:22 

0.0439±0.0031 

0.330±0.010 

1.399±0.042 

5.048±0.053 

1.14285 

0.53285 

-61.75 

Spitzer-IRAC 

2007-05-08 

19:50 

0.0492±0.0041 

0.356±0.009 

1.540±0.042 

5.203±0.045 

1.14291 

0.53295 

-61.74 

Spitzer-IRAC 

2007-05-08 

20:16 

0.0486±0.0042 

0.317±0.009 

1.370±0.044 

4.907±0.052 

1.14297 

0.53304 

-61.73 

Spitzer-IRAC 

2007-05-08 

20:39 

0.0447±0.0031 

0.316±0.009 

1.435±0.044 

5.018±0.038 

1.14302 

0.53312 

-61.73 

Spitzer-IRAC 


UT 

Wavelength 

(/rm) 

Flux 

(mJy) 

'^helio 

(AU) 

^obs 

(AU) 

a. 

(°) 

Observatory 

Instrument 

2011-09-09 19:01 

70.0 

27.0±1.7 

1.0146021 

0.1742221 

+85.3 

Herschel/PACS 

2011-09-09 19:59 

70.0 

24.3±1.5 

1.0144716 

0.1742140 

+85.4 

Herschel/PACS 

2011-09-09 19:25 

100.0 

14.2±1.1 

1.0145368 

0.1742180 

+85.4 

Herschel/PACS 

2011-09-09 20:24 

100.0 

12.0±1.2 

1.0144064 

0.1742101 

+85.4 

Herschel/PACS 

2011-09-09 19:25 

160.0 

6.5±2.3 

1.0145368 

0.1742180 

+85.4 

Herschel/PACS 

2011-09-09 19:59 

160.0 

4.4±1.9 

1.0144716 

0.1742140 

+85.4 

Herschel/PACS 

2011-09-09 20:24 

160.0 

7.6±2.0 

1.0144064 

0.1742101 

+85.4 

Herschel/PACS 

2011-09-17 09:21 

8.59 

23.7±3.6 

0.9908803 

0.1798955 

+89.4 

ESC VLT/VISIR 


radar-resolved shape model llNolan et al.ll201^ for Bennu in 
our fitting procedure, and the incorporation of that shape 
model will definitely p rovide an improved determination 
of th e thermal inertia llEmerv et alJI20l3 : IYu. Ji fc Wand 
l2014h . 

According to iFowler fc Chillemil ll 19921 ). an asteroid’s 
effective diameter Defi, which is defined by the diameter of 
a sphere with a volume identical to what the radar-derived 
shape model describes, can be related to its geometric albedo 
p„ through its absolute visual magnitude by the following 
equation: 


DeS = 


1329 X 10-^’'/® 


(km) . 


( 3 ) 


Thus, we actually have three free parameters - thermal iner¬ 
tia, roughness fraction and effective diameter (or geometric 
albedo)- that should be extensively investigated in the fit¬ 
ting process. Other parameters are listed in Table 

Herein the surface roughness is modeled by a frac¬ 
tional coverage of hemispherical craters, symbolized by 
/r, whereas the remaining fraction, 1 — /r, represents a 
smooth fiat surface on the asteroid. The hemispherical crater 


adopted in this work is a low-resolution model consisting 
of 132 facets and 73 vertices, following a t r eatment simi¬ 
lar to that shown in (i Rozitis fc GreenI [20111 : IWolters et al.1 
I2OIII : IYu. Ji fc Wandl2014h . As the sunlight is more easily 
scattered over a rough surface than a smooth flat region, 
roughness can decrease the effective Bond albedo. Using the 
above mentioned model of surface roughness, the effective 
Bond albedo Agg of a rough surface can be relevant to the 
Bond albedo A b of a smooth flat surface and the roughn ess 
fraction /r by (IWolters et ahllioil] : IYu. Ji fc Wang|[20l4) 

Aeff = /r 2 _ + (1 - /r)^b . (4) 

On the other hand, the effective Bond albedo Aes is related 
to geometric albedo by 

Aeff = Pv^iph , (5) 

where is a pha se integral that can be approximated by 
llBowell et al.lfl98^ 

Pph = 0.290 -b 0.684G , (6) 

where G is the slope parameter in the H, G magnitude sys- 
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Table 2. Assumed physical parameters used in ATPM. 


Property 

Value 

References 

Number of vertices 

1348 

fNolan et al. 2013) 

Number of facets 

2692 

fNolan et al. 2013) 

Shape (a:b:c) 

1.1135:1.0534:1 

(Nolan et al. 2013) 

Spin axis 

(-88.0° ,45.0°) 

fNolan et al. 2013) 

Spin period 

4.2976 h 

("Nolan et al. 2013) 

Absolute magnitude 

20.40 

fHererenrother et al. 2013) 

Slope parameter 

-0.08 

IHergenrother et al. 20131 

Emissivity 

0.9 

assumption 


tern of iBowell et alJ lll989h . Then for each thermal inertia 
r, roughness fraction /r, and effective diam eter Deft case, 
a flu x correction factor FCF is defined as l|Wolters et alJ 

[2^ 

FCF = ^ 

1 ,initial 

where AB,now is calculated by inversion of equation (0, to fit 
the observations, and then the so-called reduced defined 
as dMuller et alJlioilI l 


Xreduced 


n 


FCF ■ Fmadel “ -F’obs(Ai) ' 


CA, 


( 8 ) 


can be obtained to assess the fitting degree of our model 
with respect to the observations. Herein the predicted model 
flux Fniodei(r, fn, Detf, A) is a rotationally averaged profile, 
since the rotation phase of Bennu was unknown at the time 
of observation. In addition, the FCF plays a less significant 
role in determining thermal inertia, but simply brings about 
~ 0.2% influence on the outcome of effective diameter and 
geometric albedo. 

In order to simplify the best-fit searching process, 
a set of thermal inertia are given in the range 0 ~ 
1000 Jm“^s“° ®K“^. For each case, a series of roughness 
fractions and effective diameters are evaluated to find out 
which can be taken as the likely solution with respect to 
the observations. The fitting outcomes are summarized in 
Table |31 In this table, the values are thereby relevant 
to each thermal inertia, roughness fraction and effective di¬ 
ameter. Each effective diameter listed in Table 0 is a pro¬ 
file that gives a smallest for each thermal inertia and 
roughness fraction pair. Roughly speaking, the values 
imply a best-estimate solution of thermal inertia between 
200 ~ 300 Jm“^s“°'®K“^, which can be derived from the 
smallest x^ value in the two-dimension phase space of x^ 

(r, /r). 


3 RESULTS ANALYSIS 

3.1 Thermal inertia and Roughness fraction 

To obtain the best-fitting solution of thermal inertia from 
Table 0 the F ~ curves are plotted to determine how 
globally changes with the free parameters of thermal inertia, 
roughness fraction and effective diameter (see Figure [T|. 

In Figure [U the bold black curve represents the best 
estimated solution in the fitting, which is a cubic spline in¬ 
terpolation curve for each lowest y^ calculated from each 



Figure 1. F ~ y^ profile fit to the observations. Each dashed 
curve represents a roughness fraction /r in the range of 0.0 ~ 1.0. 
The bold black line is a cubic spline interpolation curve for each 
lowest y^ derived from each free parameter. 



Figure 2. y^ (F, /r) contour according to Table |3] The color 
(from blue to red) means the increase of profile of y^. The blue 
curve labele d by y^;„ -fl is assumed to be Itr limit to the free fit 
parameters llEmerv et al.ll^l4l : iBevington &: Robingsonll2003h . 


thermal inertia and roughness fraction. From the curve, we 
find Xmin (referred to the minimum y^), occurs in the case 
of F = 240 Jm“^s“° ®K“^, /r = 0.04, and Uefi = 510 m. In 
addition, we note that y^in decreases rapidly when F is less 
than 240 Jm~^s“°'®K“^, but increases rather slowly when F 
becomes larger than 240 Jm“^s“°'®K~'^. This phenomenon 
results from the absence of observations at low phase angle, 
and thus places limitations on the parameter space of ther¬ 
mal inertia. However, the combination of observations at 
several different phase angles actually weakens the above- 
mentioned disadvantage. Therefore, we can still achieve the 
best estimate solution for thermal inertia and roughness for 
Bennu. 















































Table 3. ATPM fitting results to the observations.(The combined effective diameters DefF ^re given in m) 


Roughness 

fraction 

/r 

0 

Deft 

Thermal inertia V ( Jm '^■®K 1) 

100 200 300 

^efl ^eff DeS X^ 

400 

.Deff X^ 

500 

X^ 

600 

DeS X^ 

800 

Deft 

1000 

X^ 

0.00 

475.1 

87.04 

496.2 

17.87 

516.8 

4.97 

535.2 

9.49 

550.6 

20.00 

563.3 

31.85 

573.7 

43.55 

589.9 

63.77 

601.1 

79.70 

0.05 

457.1 

98.04 

478.0 

24.37 

497.6 

6.48 

514.8 

5.85 

529.3 

11.51 

541.1 

19.06 

550.8 

26.96 

565.9 

41.08 

576.3 

52.50 

0.10 

441.0 

108.07 

461.6 

31.04 

480.2 

9.45 

496.5 

5.00 

510.0 

7.08 

521.0 

11.48 

530.1 

16.60 

544.1 

26.32 

553.8 

34.47 

0.15 

426.6 

117.23 

446.7 

37.66 

464.5 

13.29 

479.8 

5.96 

492.5 

5.36 

502.8 

7.40 

511.3 

10.46 

524.3 

16.90 

533.4 

22.64 

0.20 

413.5 

125.59 

433.2 

44.11 

450.2 

17.63 

464.6 

8.09 

476.5 

5.45 

486.2 

5.69 

494.2 

7.17 

506.3 

11.15 

514.8 

15.05 

0.25 

401.5 

133.25 

420.9 

50.30 

437.1 

22.20 

450.8 

10.97 

462.0 

6.74 

471.1 

5.58 

478.6 

5.86 

489.9 

7.93 

497.9 

10.42 

0.30 

390.6 

140.26 

409.5 

56.21 

425.1 

26.87 

438.1 

14.30 

448.7 

8.82 

457.3 

6.58 

464.3 

5.90 

474.9 

6.47 

482.4 

7.85 

0.35 

380.5 

146.72 

399.1 

61.81 

414.0 

31.52 

426.4 

17.89 

436.4 

11.43 

444.5 

8.31 

451.2 

6.87 

461.2 

6.25 

468.3 

6.73 

0.40 

371.3 

152.66 

389.5 

67.12 

403.8 

36.08 

415.6 

21.62 

425.1 

14.36 

432.8 

10.53 

439.1 

8.48 

448.6 

6.88 

455.2 

6.64 

0.45 

362.7 

158.16 

380.5 

72.14 

394.3 

40.53 

405.6 

25.39 

414.7 

17.48 

422.0 

13.07 

427.9 

10.52 

436.9 

8.13 

443.2 

7.28 

0.50 

354.8 

163.24 

372.2 

76.88 

385.5 

44.82 

396.3 

29.14 

405.0 

20.70 

411.9 

15.81 

417.6 

12.85 

426.1 

9.78 

432.1 

8.43 

0.55 

347.4 

167.97 

364.5 

81.36 

377.3 

48.95 

387.7 

32.83 

395.9 

23.95 

402.6 

18.66 

408.0 

15.35 

416.1 

11.72 

421.8 

9.94 

0.60 

340.5 

172.36 

357.2 

85.58 

369.6 

52.92 

379.6 

36.44 

387.5 

27.18 

393.9 

21.56 

399.0 

17.95 

406.8 

13.84 

412.3 

11.70 

0.65 

334.0 

176.46 

350.5 

89.57 

362.5 

56.72 

372.0 

39.94 

379.7 

30.38 

385.7 

24.47 

390.7 

20.60 

398.1 

16.08 

403.4 

13.62 

0.70 

328.0 

180.29 

344.1 

93.35 

355.7 

60.35 

365.0 

43.32 

372.3 

33.50 

378.2 

27.35 

382.9 

23.27 

390.0 

18.38 

395.0 

15.65 

0.75 

322.4 

183.87 

338.2 

96.92 

349.4 

63.82 

358.4 

46.59 

365.4 

36.55 

371.0 

30.18 

375.6 

25.91 

382.4 

20.71 

387.3 

17.73 

0.80 

317.1 

187.24 

332.6 

100.30 

343.5 

67.13 

352.2 

49.74 

359.0 

39.51 

364.4 

32.95 

368.7 

28.52 

375.3 

23.04 

380.0 

19.84 

0.85 

312.1 

190.40 

327.3 

103.50 

338.0 

70.30 

346.3 

52.77 

352.9 

42.37 

358.1 

35.65 

362.3 

31.07 

368.6 

25.34 

373.1 

21.95 

0.90 

307.4 

193.38 

322.4 

106.54 

332.7 

73.32 

340.8 

55.68 

347.2 

45.14 

352.2 

38.28 

356.3 

33.56 

362.4 

27.61 

366.7 

24.04 

0.95 

303.0 

196.18 

317.7 

109.43 

327.8 

76.21 

335.7 

58.47 

341.8 

47.81 

346.7 

40.82 

350.6 

35.99 

356.5 

29.84 

360.7 

26.11 

1.00 

298.9 

198.84 

313.4 

112.18 

323.2 

78.97 

330.8 

61.15 

336.8 

50.38 

341.5 

43.28 

345.3 

38.35 

351.0 

32.01 

355.1 

28.14 
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Table 4. Derived properties of Bennu from ATPM. 


Property 

Result 

Thermal inertia F ( 

240+^« 

Roughness fraction /r 


Effective diameter 


Geometric albedo pv 

u.U4^_0_ooil 


Figure [2] shows a contour of in the 2-dimensional pa¬ 
rameter space (r, /r), in which is represented by colour. 
The increase of is shown by ColorBar from blue to red. 
The black shows where Xmin occurs in the (F, /r) param¬ 
eter space. The blue curve corresponds to Xmin + where a 
la limit of the free fit parameters F and /r is assumed. We 
observe that the blue profile is closed in the (F, /r) space, 
indicating that we can provide a constraint for thermal in¬ 
ertia and roughness fraction simultaneously, within the la 
limit. However, when considering a higher limit of free pa¬ 
rameters above la, the degeneracy of thermal inertia and 
roughness fraction cannot be removed as well as Icr level, 
indicating that thermal inertia and roughness fraction may 
be simply separated at la level based on the calculations 
given in Table O Therefore, if la limit is reliable, we may 
safely conclude that the roughness fraction is likely to be in 
the range of 0~0.3, whereas the thermal inertia is possibly 
in the range of 180 ~ 680 Jm“^s“° ®K“^. Our result coin- 
cides with earlier i nvestigations of iMitller et al.l (l2012f l and 
lEmerv et al.l (l2014lf . Table |4] summarizes the derived results 
of the properties of Bennu. 

In the subsequent section, we will now employ these 
derived parameters to evaluate the surface thermal environ¬ 
ment of Bennu at its aphelion and perihelion, respectively. 


3.2 Temperature Distribution 

Since Bennu is a spacecraft target asteroid, the surface tem¬ 
perature is very essential to be estimated as a reference for 
the mission input. Now the surface nature of Bennu’s ther¬ 
mal inertia and geometric albedo are precisely derived from 
the thermal modelling process, therefore we can easily ob¬ 
tain the surface temperature distribution of Bennu at any 
epoch from the thermophysical model. 

Figure [3] shows the global surface temperature distri¬ 
bution of Bennu at its aphelion (left panel) and perihe¬ 
lion (right panel), respectively. In the Figure, the utilized- 
coordinate can be named as ’’the frame system”, where the 
coordinate origin represents the asteroid center, the z-axis 
coincides with the positive spin axis, and the x-axis follows 
the rules so that the Sun always locates within the x-z-plane. 
The profile of temperature in Figure [3] is shown by the col- 
orbar index - the red region represents the facets are sunlit, 
whereas the blue facets are related to relatively low temper¬ 
atures. As shown in Figure [31 the surface temperatures of 
its aphelion and perihelion are roughly in the range of 160 - 
300 K and 150 - 400 K, respectively. 

Figure 2] shows the equatorial temperature distributions 
of Bennu at its aphelion and perihelion, respectively. The 
maximum temperature does not appear at the sub-solar 


point, but delays ~ 28°, and the minimum temperature oc¬ 
curs just a little after the local sunrise, trailing ~ 12°. This 
postponed effect between absorption and emission is actu¬ 
ally brought about by non-zero thermal inertia and the finite 
rotation speed of the asteroid. On the other hand, Figure [4] 
shows that the equatorial temperature of Bennu alters from 
220 to 400 K over an entire orbital period. 


3.3 Regolith 


As mentioned previously, Bennu has been chosen as the tar¬ 
get asteroid of the OSIRIS-REx sample return mission; thus 
we show great interest in the surface features of the aster¬ 
oid, whether a regolith layer exists on its surface. Generally, 
thermal inertia is a good indicator to infer the presence or 
absence of loose material on the asteroid’s surface. As is 
well known, fine dust has a very low thermal inertia of ~ 
30 Jm~^s“° ®K“^, and lunar regolith corresponds to a rel¬ 
atively low value about 50 Jm~^s“° ®K“^. In comparison, 
the soil of a sandy regolith like Eros may have a value of 
100 — 200 Jm“^s“°'®K~^, but coarse sand (e.g., Itokawa’s 
Muses-Sea Regio) is likely to possess a relatively higher 
thermal inertia profile ~ 400 Jm~^s“'^ ®K~^. However, bare 
rock has an extremel y high thermal in ertia, greater than 
2500 Jm~^s“°'®K~^ ( Delbo et al.ll200^ F In summary, the 
given information suggests that a lower thermal inertia for 
the asteroid may be related to a regolith layer. Since the de¬ 


rived thermal inertia F = 240ltfin° Jm 


-2„-0.5t^-1 


K ^ of Bennu 


is larger than that of Eros but lower than that Itokawa, we 
may infer that regolith may exist on the surface of Bennu. 

Although thermal inertia is associated with the sur¬ 
face properties, a question arises - what does the profile 
of thermal inertia tell us about the surface properties? Gen¬ 
erally, we are interested in the investigation of grain size 
and regolith thickness of the asteroid, because these play 
an important role in understanding the primitive substance 
on/underneath the regolith layer; they further provide good 
engineering parameters for the sample return mission. The¬ 
oretically, the so-called ’skin depth’: 


h = 



pc^/uj’ 


(9) 


is usually presented to characterize the maximum grain size 
of regolith. Thus, if the k, p, c and w of Bennu are known, 
in general Is can be estimated to reveal the maxi¬ 
mum grain size of regolith on the surface of Bennu. 
Recently, the bulk density of Bennu ha s been updated to 
1.26 ± 0.07 g cm“® llGheslev et al.|[2014^ . However, the re¬ 
golith density for Bennu is unknown, so the bulk density 
may work as a reference because the regolith density is gen¬ 
erally no larger than the bulk density. Thus we adopt the 
bulk density as an approximation for the average density of 
the surface regolith when estimating the skin d epth Is- The 
rotati on period is about 4.2976 h according to lNolan et al.l 
ll2013ll . We may use the specific heat capacity of CM or Cl 
carbonaceous chondrites to approximate the specific heat 
capacity of the su rface regolith, where c ~ 500 Jkg“^K“^ 
liOneil et al.ll201^ . In this way, we can estimate the skin 
depth of Bennu to be about 1.9 cm, which suggests that the 
grain size may be less than cm-scale. Of course, the estima¬ 
tion of grain size from skin depth is rough; thus it would be 
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Figure 3. Global surface temperature distribution simulated by ATPM based on the derived thermal inertia F = 240 Jm~^s~®'^K“^ 
at the aphelion and perihelion, respectively. The colorbar indicates the range of temperature, where red for high temperature and blue 
for low temperature. 


Equatorial temperature 



longitude (°)(0° is subsolar point) 


Figure 4. Equatorial temperature distribution of Bennu at aphelion and perihelion, respectively. 
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appreciated if another way could be found to estimate the 
grain size of Bennu’s surface regolith. 

According to the definition of thermal inertia: 

r = ^(pPpCK , (10) 

where p is mean density of the surface, Pp is the grain’s den¬ 
sity, (j) (1-porosity) represents packing fraction of surface, c 
is the specific heat capacity and k is the thermal conductiv¬ 
ity. Thus the thermal conductivity k is directly related to 
the thermal inertia F and pac king fraction cj>. 

iGundlach fc Bluinl (l2013l ) builds a thermal conductivity 
model and provides a formula to theoretically estimate ther¬ 
mal conductivity k from average radii of grains r, packing 
fraction (f) and temperature T : 

K{r,T,cf>) = Ksoiidl^ ^ ^ ■ [/ie2:p(/2<(>)] 'X 

-fSo-er^ei ^ . (11) 

4 > 

The details of eq. dill) are described in (IGundlach fc BlurnI 
l2013f) . Thus, the relationship between thermal inertia T and 
grain size r can be obtained by comparing the k derived 
from r to the k estimated from Formula ini). Hence, we 
can estimate the surface grain radius r with the thermal 
inertia derived from above-mentioned thermophysical mod¬ 
elling process. 

Figure [S] shows that the r ~ F curve is plotted with 
the parameters listed in Table [5] in combination of Equa¬ 
tion HQ] and HD We assume an average surface tempera¬ 
ture T = 300 K based on Figure [4] in the computation. 
The o ther parameters are adopted from IGundlach fc Blrirnl 
ll2013li . Since there are still many assumptions and uncer¬ 
tainties in the parameters of Table [S] the uncertainties of 
(El) are rather difficult to determine. Thus we just choose 
the suggested parameters to estimate the most likely re¬ 
sults for Bennu. Using the best-fit value of thermal inertia 
F = 240 Jm“^s“° ®K“^, we estimate the grain radius is 
likely to be in the range 2 ~ 5 mm. Furthermore, consid¬ 
ering a la range of thermal inertia, we may estimate the 
grain radius is possibly in the range between 1.3 ~ 31 mm. 
According to this evaluation of grain radius, we infer that 
boulders or rocks may be few on the surface of Bennu, im¬ 
plying that the Touch-And-Go Sample Acquisition Mecha¬ 
nism (TAGSAM) designed by the OSIRIS-REx team will be 
available to be implemented successfully. 


4 DISCUSSION AND CONCLUSION 


Table 5. Physical parameters in eg. ll II (IGundlach fc Bluml2013^ . 


Property 

Value 

^solid 

1.19 -1- 2.1 X IQ-^T [Wm-iR-i] 

h 

0.25 

E 

7.8 X IQio [Pa] 

1{T) 

6.67 X 10-®r [Jm-2] 

h 

0.0518 ±0.0345 

h 

5.26 ±0.94 

X 

0.41 ± 0.02 

ei 

1.34 ±0.01 

e 

1 

Pp 

3110 [kgm-3] 

c 

560 [Jfcg-iA-i] 

T 

300 K 


Kgolid: thermal conductivity of the solid material 

p: Poisson’s ratio 

E: Young’s modulus 

7 (T): specific surface energy 

e: emissivity of the material 

Pp -. density of the solid material 

c: heat capacity of the solid material 

T: surface temperature 



Grain size (mm) 


Figure 5. r ^ T cu rve of Bennu based the method of 
iGuudlach &; BlumI l|2Q13h . giving the possible range of grain ra¬ 
dius r between 2 5 mm for the best fit thermal inertia 

240 While the la range of P is considered, the 

possible grain radius r may be between 1.3 ~ 31 mm . 


In this work, using the thermophysical mo del incorpo¬ 
ratin g the 3D radar-derived shape model (| Nolan et alJ 
l2013lf . and mid-infrared data of Spitzer-PUI, Spitzer- 
IRAC, Herschel/PACS and ESQ VLT /VISIR observations 
dMiiller et al.l[20ll : lEmerv et al.ll20l3 f , we have investigated 
the thermophysical nature of Bennu. We show that the ther¬ 
mal inertia F = 2401*1^0° Jm“^s“°'®K“^, the roughness frac¬ 
tion /r = 0.04jlQ g 4 , the effective diameter Deff = 5101 |q m, 
and the geometric albedo pv = O.OdTjlQ Qg®,. The effective 
diameter ac quired herein is a little larger than the for¬ 
mer results dMiiller et al.l[20ll : lEmerv et al.l 201^1 and that 
of the radar-derived shape model dNolan et al.l 20131 '). but 
the value of geometric albedo we derived is much closer 


to those of iMiiller et al.l d2012l l and lEmerv et all d2014|j . 
because a greater absolute visual magnitude is used in 
the fitting process. The be st solution of therm al inertia 
is even lower than that of lEmerv et al.l d2014ll . whereas 
the low roughness results are consistent with each other. 
The slight dev iation of Bennu’s t hermal inertia between 
this work and lEmerv et al.l d2014l ~) may arise from three 
aspects: first of all, we adopt a different thermophysical 
mode lling procedure dRozitis fc Greenll201ll : IYu. Ji fc Wand 
l2014h . where the beaming is modelled by solving the energy 
conservation equation in consideration of multiple scatter¬ 
ing of sunlight and self-heating of thermal emission; then. 
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the model-deri ved fluxes are rotat ionally averaged in this 
work, whereas lEmerv et alJ ll20l4l constrained the rota¬ 
tional phase directly from Spitzer data; last, the thermal in¬ 
ertia and roughness fraction are simultaneously determined 
in the 2D parameter space, thus lea ding to a larger poss ible 
range of thermal inertia than that of lEmerv et alJ (|2014l 'l. In 
a word, the physical parameters for Bennu w e deriv ed are 
essentially supportive of those of lEmerv et alJ (120141 '). 

We em p loy t he r atio of ’observation/model’ 
l|Muller et al.l 1200,1 l201ll. l2012ll to examine how the 
theoretical model results match the observations at various 
phase angles and wavelengths (see Eigure [S] and O, for the 
reliability of our fitting process and derived outcomes may 
be verified from these comparisons. 

In Figure [6l the observation/ATPM ratios are shown 
at each observational wavelength for /r = 0.04, E = 
240 Jm~^s“° ®K~^, and Defr = 510 m. The ratios are dis¬ 
tributed nearly symmetrically around 1.0, despite several 
ratios at low-wavelength 3.6 pm and long-wavelength about 
100 pm move relatively farther from unity. From Figure 
[T] we find that the deviations of low-wavelength and long- 
wavelength come from Spitzer-IRAC and Herschel/PACS, 
respectively. As the Wien peak of the thermal emission of 
NEAs approximately arises at 10 pm, the deviation at those 
wavelengths far away from the Wien peak is inevitable, 
indicating that He rschel observations at long-wavelength 
dMilller et al.l[2012l l are not very sensitive to the averaged 
surface temperature of Bennu, and thus have insignificant 
influence on the derived thermal inertia. However, the Her¬ 
schel data can still be helpful to improve the determination 
of the effective diameter. 

On the other hand, although the Spitzer-IRAC and Her¬ 
schel/PACS observations appear to deviate away from unity, 
the degeneracy of thermal inertia and roughness are actually 
removed via our modelling procedure through the combina¬ 
tion fitting to all observations (see Table [!} at discrepant 
phase angles. Therefore, we successfully provide a constraint 
for thermal inertia and roughness fraction simultaneously 
based on the Icr limit. This processing method differs from 
previous models, which usually determine thermal inertia 
with several empirical roughnesses. Hence, we safely come 
to the conclusion that the outcomes of thermal inertia and 
roughness are reliably derived from our modelling process. 

In conclusion, the average thermal inertia of Bennu of 
approximately F = 240/lgQ° Jm“^s“° ®K~^, provides direct 
evidence that the asteroid’s surface properties are an inter¬ 
mediate case between Eros and Itokawa, implying the possi¬ 
ble existence of regolith on its surface. As described above, 
a very small roughness fraction for the asteroid is derived 
simultaneously with thermal inertia in the modelling proce¬ 
dure, therefore we infer that Bennu’s surface would be fairly 
smooth. Additionally, based on the best-estimate thermal 
inertia and roughness, fine-grained regolith would be likely 
to exist and cover a large area of the surface of Bennu. Fi¬ 
nally, the estimated grain size ranging from 1.3 to 31 mm is 
indicative that Bennu would be an excellent target asteroid 
for the forthcoming OSIRIS-REx sample return mission. 



Figure 6. The observation/ATPM ratios as a function of wave¬ 
length for /r = 0.04, r = 240 and Deft = 510 m. 


2.0 I-r 




A 
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Spitzer-PUI 
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Figure 7. The observation/ATPM ratios as a function of phase 
angle for /r = 0.04, P = 240 , and Deft = 510 m. 
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